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ABSTRACT 


A statement of a method [1] for obtaining an approx- 
imate lower 100(l-a)% confidence limit on System Reliabil- 
ity is made. A discussion of evaluating the accuracy of 
this method on a digital computer is presented. Following 
this the development of a continuity correction factor is 


developed and the accuracies of this refinement for 


= 


several hypothetical systems are presented and discussed. 
eis 
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CHAPTER I 
PURPOSE OF THE STUDY 

In past years our advancing technology has created 
systems which are becoming more and more complex. In the 
military this has been particularly true, and thus, has 
produced a great need for a method of predicting the over- 
all reliability of these systems. Many problems in the 
area of reliability have not been solved. One of the 
unsolved problems is the formulation of a method for compu- 
ting a confidence interval on system reliability where the 
reliability of the system is computed as the product of the 
reliabilities of each of the components, and the sizes of 
the test samples of the components are not all equal. If 
the test samples are equal, there exist well-known methods 
for computing a confidence interval. Because of the size 
and complexity of present day systems it is necessary for 
major contractors to have many sub-contractors. They may 
be scattered througout the country and may produce a limited 
Supply of each component, of which very few may be used for 
reliability testing. These factors plus the general practice 
of testing before the system is assembled points out that 
having the sample sizes all equal is a rare event. 

An ad hoc procedure for computing an approximate con- 
fidence interval has been developed, based upon Bernoulli 
test data from unequal samples of the various components of 
the system [1]. This technique will be called Method I 


throughout this document. 


An underlying assumption to this procedure is that the 
distribution of the unbiased estimator of the negative 
natural logarithm of system reliability may be approximated 
by a two-parameter Gamma probability distribution. Conse- 
quently, a discrete set of values are approximated by a 
continuous function. Realizing the inherent error in this 
approximation this Thesis has as its purpose 
(1) to develop a correction factor to reduce this 
Srner: and 

(2) to simulate the distribution of the lower con- 
fidence limit of various systems and examine 
the results achieved with and without this 
correction. 

The fruitfulness of the investigation that follows is 
demonstrated by a discussion of the results. In later 
chapters, it will be explained that the a th percentile 
of the simulated distribution of the lower confidence limit, 
called Rs ua)’ and the true reliability of the system, Reo 
should be the same; therefore, the difference will be used 
as a measure of the accuracy of Method I. Consider the 
Following example: a system is made up of 15 components; 

14 have a sample size of 20 to be tested; and each a failure 
probability of .995; one has a sample size of 150 and prob- 
ability of failure .850. The true reliability, Ro» is..7924. 


The ao th percentile point of the distribution of the lower 


100(1l-a)% confidence limit, without correction were computed 


to be: 


4 ae) IRs ,L(a) s | 
S56 8170 0246 
10 8175 0251 
20 8000 0076 


Now by adding the correction factor the results became: 


OL Re (a) Peet io), 
a0i5 . Silas .01910 
Pel 0) . Saleas 20 2a 
20 .7969 .0045 


The above example certainly did not have the best 
accuracy of the various systems chosen nor was it the 
worst. The following summary shows the average and vari- 
ance of the error in all the cases studied. 


Without correction: 


a Average Error Variance of Error 
.05 .0438 00213 
-10 0305 .00076 
.20 0109 .00016 


With correction: 


(o} Average Error Variance of Error 
20S .0286 .00041 
- 10 .0168 ~O001L5 
~20 -0056 -00005 


The next chapter plus the appendices contains a com- 
plete explanation of Method I and discusses some further 
approximations made in the method. 

Chapter three discusses the digital computer's role in 


obtaining the distribution of the lower confidence limit of 


the reliability by Method I. Also mentioned is the measure 
of accuracy used to evaluate the method's precision. 

Chapter four outlines the reason for the correction 
factor and discusses the logic behind its development and 
application to this method. 

In Chapter five the results of various computations 
are presented and discussed. 

The last chapter summarizes the study and states the 


CONCLUSIONS, 
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CHAPTER I1 
EXPLANATION OF METHOD I [1] 
If a system consists of k components in logical series, 


the true system reliability, Ro, may be expressed as follows: 


(1) 


where Pp, is the true reliability of the aoe component. The 
component may be either continuous-operating or cycle type. 


It is desired to find a lower 100(1-a)% confidence 


interval, for R.. That 15 pea, statistic such that 
Ss s,L(a) 

whatever the actual values of Pur Poree++Pye 
me. 2 Rs (a)! can G 


The method proposed is to put n; items of component i 
on test, i=l, 2, ...., kK; under the environmental conditions 
defined in the mission and let each operate until failure 
or the mission time is reached, whichever occurs first. 
Letting 
| f. = number of components of type i 

that didn't complete their mission. 


We define: Bs = 1-q. where q; = 


(2) 


f 
ile: 
ak 


and R = 


ON 
Ss Pi 


ian 


i=l 
The random variable =1nR_ is now used to obtain the 

lower confidence bound on R.- This is accomplished by 

approximating the distribution of — Lae by a two parameter 


gamma distribution and obtaining the confidence bound via 


the gamma distribution. 
Ihde 


To continue the development define: 
k 
— -1nR. = =) in (1-q,) (3) 
1=1 


and, expanding the natural logarithm in an infinite series 


k 
Ses tag) Flea)? $-ag)?. | 
i=l 
k 00 5 
a a (4) 
i=l j=1 > 


If each qs is small, say, .15 or less, the above series 
can be approximated by the first two terms of the series, 


called T. for ease of expression. That is, 
2 


k Ce 
= = 
by as "2 i 


T, (5) 


(ian 


Op) 
Ile 


i=l 


In Appendix III [4], it is shown that the error due to the 


apove truncation is quite small; in fact; 


k q 3 


k 
S < Ta) : (6) 

Se) eee 
In Appendix I [4], an unbiased estimator for T, is devel- 


oped and is called Wee 


2 
n . 44 
Cp Chas em 7) 
where 
2n.-3 
a. = : (8) 
ii 2(n,;-1) 
aaa 
b. = nal (9) 
and as in (2) 
a 
qi ~ “A, 
wi 


ifZ 


An unbiased estimator is considered important here because 


the 1A are to be added to obtain: 


i (10) 


S = 
i 


| an 


il 
Therefore, a positive or negative bias on all of the sae 
would be undesirable. It also follows that § will be an 
unbiased estimator of S. 

An approximate value for the variance of § can be 


shown to be (see Appendices II and IV of [4]) 


B | hJ> 
jJ- 


bh. 


Var (T.) = 
1 + i 


ce) 


dl 


iar 


, k 
Varese) 
= 1 


The distribution of S is now fitted with a two paramr- 
eter gamma family. That is, it is assumed the probability 


distribution of § is given by the density function fa 


(x,r,9) defined as: 


1 Se; 
£,(x,r,0) a aera ,X?0,Fa6 ,220 
IP eS, 
= 0 elsewhere. 


Estimates of f and 6 are found by the Method of Moments 


and the fact that Z is a Xoy random variable is used to 
obtain a confidence interval for -Ink. and, thus, from this 
the confidence interval for Re can be derived. 
Since § is assumed to be gamma distributed: 
E[S] = Or (12) 
Var(S] = Q*r (13) 


But also from (10) that 
in k 
Eis] = ) vf. 


abs 


and 


. k 
eue (SS) 


Thus, solving the above equations for r and 0, the follow- 
ing expressions are obtained: 


2 


Ke 
il 


(14) 


et rae 


j 
T. 
es ee os 


el by 
- 


i=l 


k k 
a= ) — 4 T. (a5) 
ae T=] 


Thus, an estimator f, for r is taken as 
A 2 T. 
f= w) T ) =r (16) 


Now to obtain the confidence interval the following 
procedure is used. Since 28 is a Xan random variable and 
Or = -1nR., the following probability statements can be 


made: 


By 
o>) fe 
I'v 


2 
Xi-a,2r 


. 


Xl-o,2r 


~inR. < 2s - a 


2 
Xi-o,2r 


P -InR. < 2s ge 


vis 2 xv 
| 


2 
Xi-q,2r (17) 
It has been noticed that this last statement has an 


important property; namely 


Zi 

ze 

Xl-o,2r (18) 
is almost independent of r. As an example if a = .10 for 


14 


r = 6 the ratrermis — and for r = 12 the ratio is — 
Tt is observed that r can be varied by a factor of 2 and 
still have nearly the same value. This effect is even 


less for larger values of r. This suggests that the prob- 


ability distiwmiburion sot 


[2r] (19) where [2f] denotes:the smallest 
y g 


X1-o, [2¥] integer greater than’ or’ equal to 2? 


has a very small variance. Therefore, (18) could be re- 
placed by (19) in (17) and, thus, have the following state- 
ment 
1- a = P| -inr, < 2 4e) (20) 
X1-a, [2%] 


(20) can be rewritten as 


l-a = P R, > exp zo (21) 


X1-o, [2] 
That is a 100(1-a)% lower confidence limit for Re is 
~§ [2] 
2 
Xi-o, [24] 


= exp (22) 


RSL (a) 
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CHAPTER III 
SIMULATED DISTRIBUTION OF THE CONFIDENCE LIMiT 

In order to evaluate the accuracy of the procedure, 
a digital computer is used to simulate the distribution 
of the S's. The following parameters are used as inputs 
to the computer program: k; ns, i= 1,...., eefand P; = 1-q,, 
a a 
In order to have an Re to compare with the simulation the 
program first computes an Re fcom 

k 

RS TT, Ps : 

The program next Simulates values. A three digit ran- 
dom number is obtained using a uniform random number gener- 
ator subroutine. If the random number generated is greater 
than Pi, a failure is "counted" by the computer; if the 
random number is less than or equal to P;. no "count" is 
made. This random number generation is done n. times for 
each P;- Thus, the number of failures counted divided by 
the number of units of that component becomes an estimator 
for the unreliability of that component. 

16 
1.e., q; = — 
a 
This process is repeated for all dai thus giving an estima- 
tor for each qs Once these estimators are simulated the 
following arithmetic operations are performed by the com- 


puter and stored for further use: 
es (7) 


a6 


"> 
Il 


k ray 
ae (10) 


i= 


F> 
il 


“A 2 “A 
() T.) -_—s (23) 
ie 
i 
In order to complete the simulation of the distribution 


nw 


Ot ok replications are made of the above procedure. 


s,L(a) * 
This gives 500 values each of (10) and (23) which are used 


to compute 500 values of 


A a 2 
“INReo () = $t281/x}_4, (28) (24) 


for each value of a. These 500 points are then sorted, by 
Size, by a separate subroutine. The computer now picks out 
Ge & th percemerie of thisedistribucion, callaiteA. Then 
Sed) = et& is computed and printed out for comparison 
with the true system reliability, Re 

To develop a measure of accuracy for the simulation 
consider the underlying meaning of the probability state- 
ment 


P[R > R =a - Qg (25) 


Ss sn Wie 


If the probability density function of R is plotted 


s,L(a) 
on a coordinate system, and R, is taken to be a point on 

the abscissa, then 100(l-a)% of the area under the curve 

lies above and a% lies below R.- Therefore, R. should be 
equa! COmeee mee cl pCrecmertcr OL Chemaictmi DME1ON OF Re, L (a) 
which we have called exp (-A). 


Thus, |Ro - exp (-A) | (26) is a measure of the 


accuracy of the procedure under investigation. 
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CHAPTER IV 
DEVELOPMENT OF CONTINUITY CORRECTION FACTOR 

In the previous discussion § was assumed to be a gamma 
distributed random variable. Since a digital computer 
Simulation is used to generate the values of § it is real- 
ized that the actual distribution is discrete. Therefore, 
the program is estimating a gamma distribution. As in the 
case of estimating a binomial with a normal distribution 
a more accurate result can be had if the discrete values 
are slightly shifted. 

The following is the development of a correction 
factor (commonly called a continuity correction factor). 


Consider the following: 








yy. 
S = Aly 

TE 
; He yf Be 
Be me op | be) 

a 1 
where a. = 2n.-3 and Dare A 
a 1 1 1 
2(n,-1) n.-1 


The smallest difference between two values of S will 
be the smallest change in one of the T, terms. For a T; 
the smallest change is when f. 1s increased by one for that 
respective n. that is the max (n;) 
Mel, 2; 2 aes 
Let j} be that i 39 


ies max (ne) 
5 a 
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Aw 


and in we increase f. by one to obtain T.*; that is 


(cui) eee 
a es or 


* 
J J n. 2 n. 

J J 
Define the correction factor as 


~i (px - 
C, == @-* aT.) (29) 


where 5 


K]> 
| 
L' 

+ 
q 
| | 


(2n.-3) £.+£. 
Eee (30) 


2m. @(n .—2 
J : J 


=| ears =~ EP LG 
$3 Se ee ee ee 
NS, ns 2(nj-1) n 


Game 20) £,+2(n,-1) 
ee ee 


= : (35) 
2n. (n.-l) 
J J 


[£,°+(2n,-1)£,+2(n,-1)] - [(2n.-3)£ +£,7) 


Cc. = = je = eee «2 see ee). ae ee 
2n.(n.-l) 
= J 


Then 


2n.+2f£.-2 
ao. | a 


4n.(n.-lL 
5 J 


n.t+£ -l 
a). 


BA. - (n 1) ; (32) 
J J 


The corrected S is now defined to be 


S* = §+c 


1 (33) 


fi, 


It can be seen from the above calculations that S* 
will have a value even if all the T,'s turn out to be Zero. 
Keeping this in mind and looking at the equation for ¢f 


in equation (23) 


> 


k R 
of 3X Tie 877 J 
=| 


k n~ 
Ce—alen) eT. 
ili 

iL 


aE 


a 

b 

lt 

FY 
voile 


it is observed that in this case f becomes infinite when 
the denominator is zero. Obviously this doesn't make 
sense; therefore, a correction factor is added to the de- 
nominator. Because this correction factor has little in- 
tuitive basis other than keeping f£ from becoming infinite, 
it was picked quite arbitrarily. The following correction 


factor was chosen: 


i =D. 
Syne 
C. =e (34) 
a 
eels: 
and added to )} —_ giving a new fF ie. #*. This fagiex 
i=l “i 


proved to be adequate and didn't appear to bias the final 
results; at least, not in the region of accuracy desired 


in the model. 
The two corrected values are now combined in the prob- 


ability statement of equation (21) 
Ak A* 
l-a = P| R > exp {[-S [2 35 
5 2 exp | (22 V/ xy gy 221] (35) 
giving the associated lower confidence limit 


Ak Ak 
Re tia) = exp (-8 [22 Fs (2871 (36) 


which was simulated in the computer. 
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CHAPTER V 


RESULTS OF SIMULATIONS 


There were thirteen different combinations of param- 
eters used, representing as many hypothetical systems. 
Table I lists the input parameters and the results. Col- 
umn (1) is self-explanatory. Column (2) is the number of 
components in the system. Colums (3) gives the size of 
the sample tested for each respective component. Column 
(4) is the probability of success for each component of 
the system. Column (5) is the system reliability as com- 
puted directly from the input parameters (2), (3) and (4). 
Column (6) lists the three values of a for which a lower 
confidence limit was computed. Column (7) is the exponen- 
eral pewer of the oth percentalo-of ther drstembutvon: 

The first four cases in Table I are examples of systems 
with different component sample sizes. In Cases 5 through 
13, the sample size is constant in each case. These last 
cases show the effect that system size, sample size and 
probabilities have on the accuracy of the method. The 
effects noted will be discussed after Table II is presented. 

A further check was made on the method by computing 
Case 10 four times, each time starting the random number 
generator at a different point. The results of this test 
showed there wasnt any difference in the first four decimal 
places of the confidence limit estimate. From this, it is 
concluded the 500 replications of each case is enough to 


smooth out any fluctuations in the generated numbers. 
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TABLE es 


RESULTS OF COMPUTER SIMULATION WITHOUT 
CONTINUITY CORRECTION FACTOR 


sz | ; 3 Re oe 
ASE k i Pj Rg exp (-A) 


ee ee 0.0 




















£985 .980 
£979 .995 
1 .988 .970 .7233 
6982 .995 
.980 .968 
.967 
2 
3 
4 
5 
.7985 
7 .8325 
.8506 





Ze 


TABLE I (Continued) 


RESULTS OF COMPUTER SIMULATION WITHOUT 
CONTINUITY CORRECTION FACTOR 


(1 (2 : S, 
pa a om 





2407 
#397 + .6647 
AVA) - FOS9 





05 . 7636 
. Wile 10 , 7e52 
20 7006 


05 -6862 
~7924 - 10 . 7048 
« 210 - 7/666 








23 


Table II lists the results of the cases with and with- 
out the correction factor added. Columns (1), (2) and (3) 
were described before. Column (4) is the exponential 
pewer of thema th percentile of the distribution with ¢he 
correction added. Column (5) shows the value of the meas- 
ure of accuracy used in evaluating Method I with correction 
added. Column (6) is the same as (5) but without the 
correction added. 

From Table II, it is quite apparent, on the whole, 
that the error was reduced by the addition of the correction 
factor. It also appears the correction has less effect on 
cases where the sample sizes (n; ) vary [Cases: l, 2, 3, 4], 
but still, obviously, worth the effort. Cases 5, 6 and 7 
seemed to indicate the number of components (k) in the 
system has no effect on accuracy. Cases 10, ll, 12 and 13 
appear to point out that as the sample size (n.) increases 

wAlhG0 MC CESES 
the accuracy versus probabilities, This concludes the 
discussion of the most observable facts that can be viewed 


from the tables. 


eo. ae | 
> DC J tee words addled hy Knees. 
vot by the ea “- 


24 


TABLE II 


COMPARISON OF METHOD I WITH (*) AND 


WITHOUT CONTINUITY CORRECTION FACTOR 


SOM a hs cee ee 
“ 






.0120 .0276 | 
.0033 : .0095 
.0567 1630 | 
0227 : (0750 | 
.0177 \ Ogee aus 
1 
05 65 .0103 .0102 
; 9 .7738 a0 Teo 0086 | 0085 | 
| . 7696 .0042 , 0042 
.7495 .0429 .1062 
.7591 .0330 : .0876 
.7684 .0240 .0258 





NE OE EIRL TO tent “Ec, a 


2 


TABLE II (Continued) 
COMPARISON OF METHOD I WITH (*) AND 


WITHOUT CONTINUITY CORRECTION FACTOR 


2 3 é 5 6 
Rg QO expe (A )ameoe (-A*)-R. exp (A=A) —Re 


.05 ne ae! .0241 
11 .7924 se .7820 | .0099 
20 .7931 ‘oom 


OID 
ie2 ~ 17924 aa 
7410 


a 
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CHAPTER VI 
SUMMARY AND CONCLUSIONS 

In the preceding four chapters a procedure [1] has 
been explained which will approximate a lower 100(1l-a)% 
confidence interval on system reliability. The simulation 
of this method on a computer has been explained and the 
technique used to measure the accuracy of the method was 
stated. A continuity correction factor has been developed 
which enables the discrete computations of a digital com- 
puter to more closely approximate a continuous distribution. 

The simulation was run for thirteen cases with and 
without the continuity correction factor. These were not 
an adequate amount of cases on which to base any concrete 
conclusions. It is felt these cases are sufficient to 
show the merit of the correction factor. The following 
summary as stated in Chapter I demonstrates the validity 
of the procedure. 


Without Correction: 


) Average Error Variance of Error 
-05 .0438 .00213 
-10 0305 .00076 
-20 .0109 -00016 


With Correction: 


O Average Error Variance of Error 
. 0S .0286 .00041 
-10 .0168 00015 
.20 .0056 .00005 


From Tables I and II in Chapter V, it can be seen 


the average error for those cases where sample size (n;) 


24 


1s varied would be less than shown above. It must then 
be concluded that Method I with the correction factor 
will be an effective method for approximating the lower 
100(1 - a)% confidence limit for system reliability. 

It should be pointed out that another study; ob- 
viously, could be done in this same area. By varying the 
input parameters (k, nes p,) and determining their 
effects, one could possibly establish limits on their 
variation to keep within a desired degree of accuracy. 

In closing, I would like to express my appreciation 
to Dr. W. Max Woods of the Operations Analysis Department 
for pointing out this problem and serving as a very 
willing advisor. I also wish to extend my thanks to Dr. 
Rex H. Shudde for supplying the basic computer program 
for Method I which required slight alteration to meet 


my needs. 
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APPENDIX I 


DERIVATION OF Ty, THE UNBIASED ESTIMATOR OF tT; 


In order that 7 be an unbiased estimator of To, 
it is necessary that 


T, = E[T, ] 


Rewriting the above expression, using (5) and (7) 




















qs” aif,  b,f,° 
ena = | eee 2 
ae 
and, substituting (2) 
a," aif;  b,f,° 
qi =e uly. 
1 2n. 
1 
a 2G 2 
= — E[f.] + Eli. ] 
ne 1 2 1 
1 on. 
1 
ancl. lo 
» el a 1 2 2 
-—~n * 7. 0 4, 9, ~ 9|97> ee 
1 2n; 
«2 
s ivi a 
> “4 eee ny eke 
bd, . bd. q32 
= (a on a i on ee 
ae iL 
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Then equating coefficients of q; and —— 


al 
af 
and solving simultaneously 
ea 
a. oF i © iL. 
Ef 
om 
a gle =) 7 wae 
py 3) "3" oe Geen 
ai 
b. 
_s i. 3 
>. a om > 
- ~ 1/3 _ 
=> a, = 5(3 - b,) 
Substituting back 
Jey 
1 1. 
Dae Ca Se 
Di 
3 -b. + —= 2 
a n. 
a 
on 
i sl ~ nh. - - 
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and again substituting 


ae _ if 
ia x (3 Sa = 1? 
1 
2n. - 3 
~ 1 
~ 2(n. - l 
1 
Thus, 
2n; - 3 
T. me CW qi ag n 
1 a 


is an unbiased estimator of T.. 
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APPENDIX II 
COMPUTATION OF VARIANCE OF § 
k 


= J var (7.) 
i=l = 


“Aw 


k 
Var (S) = Var ) T 
i=l 








Weng (T.) = eee — oe [T,] 








A 2 2 
Pe 
ores 
Dos Fe el Z 
E(?.“] = El(a,q, + => )°1 
2d 
me Oe 
" 2% 2 me a re 
Ela 4 +a a; + 7 ] 
2 2 
a. a.b. lore 
_ a Z lage 3 uf 4 
i! es 4n. 
af a af 
From [3] 
2. Ze 
Elf, , sat i ay us Se er 


- 7 2) Zz," 3 
2 El(f, Eff,]) ‘lle 3E(f JEIf, ] 2E [f£,] 


Be 2.8 3 2 


- 7 2 
=n. q. 3n; qe + 2n.q; + 3niqs 


2 
= ona; ae 


™ z 4 oe 2 2 
Elf. |= EL(f£, E{f,]) ) + 4EIf, JEIf, ] 6E [f£, JE(f ]) + 
4 
3k [f.] 
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Therefore 
Z 
(2nin— 3) 
3A, 1 Ze 
Bg as mg ee ee” URE 
2° (nz =—1) n; 
(2n 2 = 3) ne 
wan temngmng {n,-q,* = 3n,“q,” + 3n,“q,” + 
ie al i i 
2(n. - 1l)n. 
at a 
Ze. 2 
3n,°q,° - 3n;q." + nq, } + 
n 
ii 4 4 3. 4 3.68 2 
ae ins q, - 6n,-q,° + 6n,°q,~ + lin,“q, 
Aine ye ne 
iL 
2 2. 2 4 3 2 
18n, q, + mm; ° ie as 6n.qs a l2n.q, mq, zt. 


N55} 


E(?;“] - 3 
ny 4 (n ) =a dik ) n Ld 
al 1 


2. 3 3 2 2 2 
(n;“q;” - 3n,q,”° + 2q,” + 3n,q," - 3q,° + q,) 


4(n, - 1) n. 
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2 2 2 
(4n. _ lan, + 9) (nq. = ao q,) + (4n. ~ 6) 


+ 


3. 4 ons 2 3 Z 3 
(nase én, 2 aloes én. Cea lin.q, 18n.q; + 


2 
4(n, - l) n. 


4 
4q,° + 4q,) + n,(1lq," + 8q,° - 24,° - 8q,) + 


> > > 
4n,(n; - 1) 


ate 2q,° a 4q;) 
2 


(-6q, * 


4n,(n; - 1) 


Let 


i 2 
B = -6q = 12q, - 4q; + 4g 
C = llq : + 8q = 2g) - 8q. 
BE i 1 u 


D = -6q,° “° 2q; f 4q; 


ele 


Then 


anos + aa + Cn. + D 
2 1 1 al 
Eee) = SS? 
ba 4n. - 8n. + 4n. 
i i i 
a, B+ 2a, (C+ 2B + 3A)n; + (D - B= 2a) 
= = FS 8 
[ 2 4n.(n. = 1) 
soe ak 
4 3 2 4 3 2 
eee on ee 
cee 0 ee 5 
al 4n.(n. - lL) 
gS 
4 | 2 
-(2q, 4q.” + 2q, ) 
a hr 
2 Z 3 2 2 
jp? fi, Si , Se es ‘Sa 
a oe 4 - + 
i n. 2n. Nn. n 2n.(n. - L) 
a i 3h si ; i 
Z 
A 2 qT q. 2 
Bic eer ee Varich ee ee eee 2, it 
nt i n. n. 
i i 
2 
2(n. - 1) 
i 


which, when simplified, gives 


2 Zs 


AD q. 3q,° = n, (2a, + 2q; - lh) 





ee Ji 
Varies + 5 — |-——— a 
1 L 7 
T; qd; 3q. 
= —_ i San Sr 1 + aoe yee - 2q. 
y Dy qu : n. i i 
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which we approximate by 


= Ty 
Var (T, ) =' a 
of 
Thus 
k T. 
Var (8) =} ae 
wy al 





=n - 2q, 
a. 


This error is further discussed in Appendix IV. 


Se 


APPENDIX III 


ANALYSIS OF THE ERROR DUE TO TRUNCATION 


tn thes transitions £rom (4) to (5) an infinite 
series was truncated to its first two terms. Examine 


the inner summation of (4) 





aa q.7 q. oa q.J 
) <a = qj a5 = ata ) ee (37) 
j= lime j=3 J 
= qT; + R, 
where R, is the remainder of the series after two terms 
; th st 
have been written out. If we test the n and (n + 1) 
terms of this series by the ratio test, 
n+l nq. 
etic Ve me ee 1 7 a 
au = 
a n+l 
n V a 
Ti 
we see that this expression is'‘always less than q. sl, 
which implies the series 1s convergent. Then bv ‘Theorer 


eis 22c7ot (2) 





ete 

1 1 

eb fil — ; 
noe" n+l Ga 


which implies that 


= n+l _ < 
oma) 7 (on + 1) e) pe 2 a. 
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Since r é qs let us choose r = qs S$ 1 so as to make 


R, as small as possible 


n+l 


R, © Gy 


n 


yf Mitte Ciencia? 


and evaluating at n = 2 
R, <q," / 30 - a (38) 
2 1 1 : 


Since q; << l in most cases of interest, choose q; = 
.15, for example, which is about the largest value it 
will ever take on, and we see that R. { 1,4 & Titan 
Substituting equation (38) into (37), 

a,” 3 

= J ee (39) 


lows 8 


j=l 


and substituting (39) into (4) 


Sas 


3 


pray 


Jl 


which upon distributing the summation sign becomes 
equation (6). For ois .05, still quite a large value 


for qs to take on, and k = 15 


< -5 
Se ae ae a oe oe LO 
ai ab 
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so the approximation 





APPENDIX IV 


INVESTIGATION OF ERROR DUDE TO APPROXIMATION IN Var (S) 


In Appendix II the Variance of the random variable § is 
computed, and it is seen that the value used in (11) before 


differs from the true variance by a factor 


2 2 
oF Eye 


i i 
2(n, =I) os n. ct ie (40) 





Ira 


i=l 
Upon embarking upon an investigation of the size of this term 
it is quickly seen that it is a task which, due to the amount 
of work involved, is beyond the scope of this thesis. A few 

general comments on this term are in order, however. 

Being a function of both the qs and ns, equation (40) 
may, for certain combinations of these parameters, actually 
be less than zero and reduce the variance of §. In order to 
determine when this occurs, it would be necessary to use a 
digital computer to compute and plot for each set of nes 
m—1,..., K, a Curve over the range of G < oes 1; and, 
Similarly, to compute and plot for each set of qa l=. oe KK; 
a curve over the range of interest of ns, 2 < n <N, where 
N might be 150. 

If the summand of (40) is plotted as a function of a 
single qa with the n. held fixed, the curve is seen to be 
Be l-Sshaped over the range 0 <q <)ly but within the real= 
istic range that q; may be expected to take on 0 < q; s ao, 


(40) is a non-decreasing function, and strictly increasing 


unless a point of inflection occurs. 


Al 


If the summand of (40) is plotted as a function of a 
Single ns, where n. is of course interger-valued, it is seen 
to be a monotone-decreasing step-function, asymptotic to the 
n. axis. 

és 
As a crude bound, on the size of the error, consider 


the summand of (40) 











2 
2 
qi 345 2 
2(n,-1) 1+ n; 7 aon 7 28l9 
2 
2 
q: | 3q. 
2 1 2 
: on. 1+ ae Aah a ei 
1 1 
2 
2 
ele 3q. 
< — a is I 
2n. ae 
1 1 
< — ‘i 
2n. n 
i i 
which implies that 
; a3" i — 2 2 ; ee 
= = - £q. ~- 4q. . >, ta — 
ee 2(n, 1) al a bea) an; n; 


which is obviously not a least upper bound, but an example 


is informative. 
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Let nu=—~20e foresal lias, ancien 


which admittedly is of an undesirable magnitude for an omit- 


ted term, but (40) is undoubtedly much less than .4. 
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